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'■ Abstract 

O 

l o j . 

The scaling limit of the less than half filled attractive Hubbard chain is stud- 
ied. This is a continuum limit in which the particle number per lattice site, 
"^J- ■ n, is kept finite (0 < n < 1) while adjusting the interaction and bandwidth 

in a such way that there is a finite mass gap. We construct this limit both 
^r) . for the spectrum and the secular equations describing the excitations. We 

find, that similarly to the half filled case, the limiting model has a massive 

00 

I and a massless sector. The structure of the massive sector is closely analo- 

gous to that of the half filled band and consequently to the chiral invariant 
SU(2) Gross-Neveu (CGN) model. The structure of the massless sector differs 
from that of the half filled band case: the excitations are of particle and hole 
type, however they are not uniquely defined. The energy and the momentum 
q I of this sector exhibits a tower structure corresponding to a conformal field 

O ■ theory with c = 1 and SU(2)xSU(2) symmetry. The energy-momentum spec- 

trum and the zero temperature free energy of the states with finite density 
coincides with that of the half filled case supporting the identification of the 
limiting model with the SU(2) symmetric CGN theory. 
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I. INTRODUCTION 



The one dimensional (ID) Hubbard model is described by the Hamiltonian 

N N , -i x , i 

H = -tJ2 ( C t,v C i+l,a + h - c ) + U J2 

i=l o-=t,J. i=l 
N 



i=l 



in which cf a (c ia ) creates (destroys) an electron at the site i with spin a, n^ a = cf a c ia . 
We impose periodic boundary conditions so the site % = N + 1 is the equivalent to the site 
i = 1. The hopping t is positive while the interaction U, for the considered attractive case, 
is negative, and the value of the chemical potential, /i, is choosen to fix the particle number 
per site, n. This model being completely integrable U and exactly solvable by the Bethe 
Ansatz (BA) plays a central role in the theory of strongly correlated electron systems 
. At the same time by direct linearization around the Fermi-points [|[ both the half-filled 
(n = 1) and the non-half-filled (n < 1) Hubbard chains can be related to relativistic field 
theory models HH, specifically to the SU(2) symmetric chiral Gross-Neveu (CGN) model. 
Since the Hamiltonian of the CGN model has been diagonalized in the continuum (with a 
somewhat unorthodox cutoff i.e. filling up the Dirac sea up to a certain depth) by the BA 
method |7|, constructing the relativistic (scaling) limit of the BA solution of the Hubbard 
chain can be of twofold interest. On the one hand it makes possible to study the relation 
between the two models in the detail, and on the other hand it promotes the Hubbard chain 
to an integrable lattice regularization of the CGN model. 

Earlier Filev ||, later Melzer ||, recently Woynarovich and Forgacs [l(J have studied 



the relativistic limit of the half-filled (HF) Hubbard chain. In ref. [ID| we have given a 



rather complete study of the structure of the states and the spectrum of the limiting model 
and also presented some convincing arguments that the scaling limit of the HF attractive 
Hubbard chain is the SU(2) (in fact SO(4)) symmetric CGN model. The scaling limit of the 
non-half-filled (NHF) attractive Hubbard chain has recently been studied by Woynarovich 
fTT| . He has shown, that the limiting model possesses a massive sector whose spectrum and 
phaseshifts agree with those found in the scaling limit of the HF Hubbard chain. 



The aim of the present work is to extend the study of Ref. [If] to the complete spectrum 
and to collect some other arguments supporting the equivalence of the scaling limit of the 
non-half-filled and that of the half-filled Hubbard chain. This also implies that the scaling 
limit of the non-half-filled Hubbard chain is the SO(4) symmetric CGN model. 

The scaling limit is a continuum limit N — > oo, a — > so that Na = L = const, (a and L 
being the lattice constant and the chain length, respectively), in which the particle number 
(N e ) per site N e /N — ► n finite, i.e. N e — > oo too. To avoid divergences, the interaction u 
has to be tuned in a special way (actually u — > 0). Finally, although a does not appear 
explicitly in H, since t oc 1/distance, t —>■ oo as a — > 0. All this is to be performed so that 



the gap in the spectrum of the unbound electrons is kept finite, that is [pf 



u — ► 0, t — > oo; N, N e — > oo, at N e /N — >• n = const < 1; 

a — >• at Na = L = const (1.2a) 
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so that 



8^ / T/ TT f 7rsin(7rn/2) 1 

mo = — y it sin {im/2) exp < > = const (1.2b) 

and 

2atsin(vm/2) = 1 , (1.2c) 

We have performed this limit both in the spectrum and in the higher level Bethe Ansatz 
(HLBA) equations of the NHF Hubbard chain. We can summarise the properties of the 
limiting theory as follows: 

(i) Like in the HF case, there are two kinds of excitations: massive and massless ones. 
While the massive sector is described in terms of well defined particles, in the mass- 
less sector the definition of the particless is not unique. (This is a major difference 
as compared to the HF case, where the excitations in both sectors are well defined 
particles.) 

(ii) The massive particles have spin 1/2. Their contribution to the energy and momentum 
is given as 

E e 0%')> and Ep( k j')> ( L3a ) 



where 



e(«) = m cosh(ft), p(n) = mo smh(n), (1.3b) 



and the rapidities. The k^s and the set of variables Xa describing the spin 

state of the particles, satisfy the following BA type equations 

Lp( Kj ) = 2*4 - J> (^^) + E ^an" 1 {^j^j , (1-4*) 



E ^an- 1 l^j^j = 2nJ a + £ 2tan~ 1 (^^) , (1.4b) 

with 

<j>{x) = -ln-^ H-) H . (1.5) 

1 r(| + if)r(i-i|) 

( |1.4| a,b) yield the same phaseshifts as those obtained for the half-filled band 

?p tr = n + (—) , ^ S = f — ) - 2 tan" 1 f — ) (Ak = Kj- - Kj-/). (1.6) 



7T / V 7T / \ 7T 
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(iii) The massless excitations correspond to states above resp. below a certain Fermi niveau 
(particles resp. holes). As the choice of the Fermi niveau is not unique, however, the 
parametrization of the states (number of particles and holes, and their momenta) is 
not uniquely defined. The energy and momentum, unlike in the HF case, consist 
not only of the sums of contributions of the individual excitations, but also contain 
certain 'collective' terms too. The spectrum of this sector shows the tower structure 
characteristic of a conformal field theory (CFT): 

E — E m = -— + — (x nm + + v ) , 
oL L 

P~Pm = Y-(s n , m + V + -U-), (1.7) 

where E m and P m being the energy and momentum of the vacuum plus the massive 
sector, are integers and 



x i 



i(n 2 + m 2 ), s ntm = nm. (1.8) 



Here the numbers n and m are integers or half integers depending on the parity of the 
number of massive particles. In the case of an empty massive sector the apexes of the 
towers yield the conformal weights (the notations are those of 



I i 

A = -(n + m) 2 , A = -(n-m) 2 . (1.9) 

coinciding with those of a c = 1 CFT with an (enhanced) SU(2) symmetry. 

(iv) The terms corresponding to the apexes of the towers in the spectrum of the massless 
sector depend nonlinearly on the number of the various excitations, hence we refer 
to them as collective terms. The parametrization of the massless sector (the Fermi 
niveau), however, can be choosen so, that the number of massive particles formally 
disappears from the collective terms and then the two sectors practically decouple. 

(v) The energy-momentum spectrum of the limiting model is the same as that of the HF 
band case: the contribution of the massless sector as a function of the parameters m 



and n coincides with the corresponding contribution found in Ref. [10], and so does the 
contribution of the massive sector too, as the parity prescription for the Ij quantum 
numbers (i.e. the quantization of the momenta) expressed in terms of the m and n 
quantum numbers is also the same as that found in Ref. [[IT]]. 

(vi) The ground state energy of states with a finite density of excitations is also in agree- 
ment with the result for the HF case. If we choose the Fermi surface so, that the two 
sectors decouple, the zero temperature free energies (which are nothing but the ground 
state energies in the presence of a chemical potential) of the different sectors become 
independent of each other, and the total free energy /(//, v) can be written as 

f(n,v) = -^n JL )-f, (i.io) 

7T mo 47T 
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where \i and v denote the chemical potentials for the massive resp. massless particles. 
The \l/(/i/mo) can be given as an asymptotic series of l/[ln(/i/mo)] for high particle 
densities. (We remark that by parametrizing the massless sector in a different way 
other terms depending nonlinearly on both chemical potentials would appear in Eq. 

The equivalence of the limiting models of the HF and NHF cases is not trivial at all. 
While in a finite length chain we see a smooth behaviour as the bandfilling n is changed, in 
the thermodynamic limit the HF and NHF cases separate. There are three major (although 
not independent) differences. 

- While the HF chain has two SU(2) symmetries (spin and isospin, this latter being 
connected to the charge), the NHF has only one. 

- The structures of the excitation spectra are different, as in the HF case the gap less 
excitations are isospin 1/2 particles, while in the NHF case these excitations are particle 
and hole type. 

- In the (naive) continuum limit of the HF case there are umklapp processes violating 
chiral symmetry, while no such terms are present in the NHF case. 

It is somewhat surprising, that some of these differences thought to be significant, disappear 



in the scaling limit: as we discussed in |10| there are indications, that the amplitude of 



the umklapp processes for the NHF chain scale out in a renormalization process, and as 
the conformal weights show, the SU(2) symmetry of the massless sector of the NHF chain 
developes in the scaling limit. The only significant difference we see after the scaling limit 
is in the structure of the massless sector. On the other hand, the fact, that the energy- 
momentum spectra of the two limiting theories coincide is a strong evidence supporting the 
equivalence of the two theories. 

The paper is organized as follows. In Sec. [H] we give the BA equations of the Hubbard 
chain, and in Sec. [Ill] we describe the solutions corresponding to the less than half filled 
band. The scaling limit of the secular equations and the spectrum are constructed in Sees. 
P"V| and |V] respectively, and the results are discussed in Sec. [VI] together with a brief review 
of the states with a finite density of excitations, and we list the differences between the 
limiting models of the HF and NHF chains in Sec. |V11| . Appendix [A] we construct the naive 
continuum limit of the model. 



II. THE BA EQUATIONS 

The eigenvalue equations of the Hamiltonian ( |1.1| ) have been reduced to a set of nonlinear 
equations by Lieb and Wu 0: 

Nk* = 2ttL - y 2 tan" 1 sinfc J ~ A ° ( 2 . la) 

M 

J, = -(modl) 
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N, 



E i A/v — sin k-j _ x— -\ i A™ — A« 

2tan~ 1 — 3 - = 2vr J a + ^ 2 tan" 1 



i=i 



Z7/4 



0=1 

A" e + M + 1 



U/2 
(mod 1) ) . 



(2.1b) 



Here N e is the number of electrons, M is the number of down spins, i.e. S z = (N e /2 — M), 
and the Ij and J a quantum numbers are integers or half-odd-integers depending on the 
parities of N e and M, as indicated. Once these equations are solved the wave-function can 
be given |T2| and also the energy and the momentum of the corresponding state can be 
calculated: 



E = NU/A - £(2t cos kj + U/2-fi), P = £ kj 
i=i i=i 



(2.2) 



For the considered U < attractive chain near the ground-state most of the electrons form 
bound pairs with wavenumbers given (up to corrections exponentially small in N) as 



sin k = A ± iu 



(2.3) 



with u = \U\/4t and A being a subset of the set A. By this relation k^ can be eliminated from 
Eq.( [2.1b| ) and one finds that the wavenumbers of the unbound electrons, the As connected 



with their spin distribution and the As of the bound pairs satisfy the equations P~5| , |T4 



2nlj = Nkj - 

n(A) n(A) 

- V?i(sin kj - X a ) - ^2 ^i(sin kj - A,,) , 

a=l T7=l 

n(A) +n(A) 



(mod 1 



(2.4a) 



n(k) n M 

(fi(K - sin kj) = 2nJ a + ^ ^{K - \p) 



3=1 



(3=1 



J ft 



n(k) +n(A) + 1 



mod 1 



(2.4b) 



with 



2-ixJ^ = N (sin^A,, - iu) + siiT^A,, + 



J n 



III 



n(k) n(A) 

<Pi(Ki - sin%) - ViiKi ~ K) 

3=1 v=l 

n(k) + n(A) + l. 
mod 1 



(2.5) 



<p m (£) = 2 tan 1 (^/um) 



(2.6) 
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Here n(k), n(A) and n(A) are the number of unbound electrons, the number of unbound 
electrons with down spins and the number of bound pairs, respectively (N e = n(k) + 2n(A), 
M = n(X) + n(A)), and the quantum numbers Ij, J a and J v are integers or half- integers as 
indicated. (As in the following manipulations the quantum numbers may be redefined by 
absorbing certain constants into them, we give the 'parity prescriptions' for the quantum 
numbers together with the equations. All these prescriptions hold for even N) The energy 
and momentum expressed by these variables is 

E = -Nut - (2i(cos kj - u) - - 
j 

- E ( 2t ( \l l - ( A ^ - m ) 2 + V 1 - + m ) 2 - 2u ) - 2 ^) > ( 2 - 7 ) 

P =^(E^-E^ + E^] • (2.8) 

\ j a V J 



III. STATES OF THE NON-HALF-FILLED HUBBARD CHAIN 

1. Ground state - the physical vacuum 

As a reference state consider the state, consisting of N r = n r (A) < N/2 bound pairs 
(and no fc's, i.e. N e = 2N r ), and the quantum number set J v is given by 

N r — 1 N r — 1 

'h — , Jrj+i — Jrj + l, JN r — ~ • (3.1) 

If both N, N r — > oo (at N r /N = n/2 kept constant) the A^'s will be distributed in an interval 
—B < A < B with a density o~ r that is defined so, that the number of A^s within the interval 
(A, A + dA) is given as Na r (A)dA. It is not hard to see from (|2.5|) , that this density is the 
derivative of the so called counting function: 

•K*)-^. 0.2) 

where 

z r (A) = — | (sin-^A - iu) + mn-\A + iu)) - — £ <p 2 (A - A„) j . (3.3) 
This relation, applying the leading term of the Euler-Maclaurin type summation formula 

hp (*) - T"f W * " ^ (/' (^) " /' (^)) («) 

r ' =1 Jx/N-XpN V V 7 V 77 

leads to the equation 
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1 f 

a r (A) = ctq(A) - — / K 2 (A - A'K(A') , (3.5) 

-B 



1 



-1 N 



a (A) = — 2R^(yi-(A-m) 2 J J , (3.6) 

K ™(0 = 1 2 Z\,2 > ( 3 - 7 ) 
{muj z + £ 2 

where the limits are determined through 

z r (±5) = ±N r /2N (3.8a) 

(i.e. 

oo — B 

2 V N 

B -oo 

as it can be seen by integrating (|3.5|) ). 

The energy of this reference state is (using summation formula fl3.4p in leading order) 



a r (A) = I a r (A) = Ul - = LJL , ( 3 .8b) 



—B 

with 



E r = -utN + J2^o( A n) 

+B 

= -utN + N J e (A)a r (A) , (3.9) 



£ (A) = - \ 4t (ibey/l - (A - iu) 2 - u^j - 2/iJ . (3.10) 

In (|3.9| ) the first term is irrelevant. The second term, iterating a r , can be transformed into 
the form 

+B 

E r + utN = N J e r (A)a (A) , (3.11) 

-B 

where e r satisfies the equation 

1 ? 

e r (A) = e (A) - — / K 2 {A - A )e r (A ) . (3.12) 

-B 

It is easy to check, that the energy of ( |3.11|) as a function of B, 1.6. cLS cL function of n, 
(at fixed fi) is minimal if 

£ r (B) = 0. (3.13) 

As £"o, depends on /i linearly so does e r , hence it is possible to choose a fi such that ( |3 . 1 3| ) 
holds at the desired bandfilling n. In the following we suppose that [i has this value, i.e. the 
reference state is the ground state. 



S 



In general the ( |3.5| ) can not be solved analytically, and B as a function of n can not 
be given in a closed form. In the relativistic limit however, when u — > 0, (|3.5| ) can be 
solved using the method devised by Yang and Yang [[L6|: the equations of the type ( |3.5| ) can 
be transformed into a series of Wiener-Hopf equations which can be solved with sufficient 



accuracy yielding [14 



Tin u I 
B = sin 1 + In 

2 7T \ 



n cos 2 ?y sin ^ ' 



2u 



(3.14) 



2. The possible excitations 

The ground state can be excited by combinations of the following three "elementary" 
excitations: 

i) introducing holes and particles in the A distribution by removing some J v from the 
ground-state set and introducing some outside it, 

ii) introducing complex As, 

hi) introducing unbound electrons (real ks and As). 

The excitations of type i) with the proper choice of fi (described above) have a dispersion 
with no gap, while those of type ii) and iii) possess gaps. The scaling limit (|1.2| ) is constructed 
so, that the gap in the spectrum of the excitations of type iii) is a fixed finite value (mo). In 
this limit the gap of the excitations of type ii) diverges, i.e. they can be discarded. (In the 
following we shall refer to type i) excitations as particles and holes or massless excitations, 
and type iii) excitations will be called massive particles or unbound electrons.) 

3. Particles and holes in the A set 

First we consider the equations of the As ( |2.5| ). To describe the particle and hole type 
excitations properly we should define the 'Fermi sea'. If n(k) + n(A) = n r (A) (mod 2) the 
Fermi sea (of Js) is the ground state set of the Js, i.e. the Fermi sea is made up by all the 
Jf integers or half-integers satisfying 

Ji = -Jf + 1/2, J f+1 = J f + l, J f <J F - 1/2, (3.15) 

with 



J F = N r /2. (3.16) 

If n{k) + n(A) ^ n r (A) (mod 2), there is no unique choice: as the parity of the J v s (i.e. that 
of the Jfs) is changed, keeping ( |3.15| ) we may chose any of Jp = (N r ± l)/2. For the sake 
of definiteness we choose 

J F = (iV r + l)/2, (3.17) 

thus 
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T Nr + 6 A _ / if n(A;)+n(A) = n r (A)(mod2) 

2 ' \ 1 if n(k)+n(A) ^ n r (A) (mod2) ' 1 j 

(We note that the Fermi sea defined above resembles most the ground-state distribution of 
the As, nevertheless many other definitions are possible: the only role of Jp is to provide a 
reference point in the space of the J quantum numbers connected to the A-rapidities, i.e. S 
can be choosen even to depend explicitly on n(k). We shall discuss this later.) 
To proceed, we define the counting function: 

z(A) = — \ (sin _1 (A„ - iu) + suT^A + iu)) - 

- - £ <Pi (A - sin kj) - - £ <^(A - A,) | , (3.19) 
in terms of which Eq. (|2.5|) reads 



z(A v ) = ^ . (3.20) 

It is clear, that ( [3.20] ) has a solution for any z(— oo) < J < z(+oo) replacing J v . The 
rapidities of the particles, holes and respectively the elements of the Fermi sea are determined 
by the equations 

z(Ap) = |, *(Aa) = ^, z(A f ) = ^, (3.21) 



respectively, where the J p s are those J v s, which are not elements of the set J/ of fl3.15|) 



while the J^s are those J/s, which are not elements of the J v set. (Note, that ±Jp themselvs 
can be neither holes nor particles, as their parity is different from that of the J v s and J/S.) 
Comparing the parity prescription J v = (n(k) + n(A) + l)/2 and the number of elements in 
( P-15| ) we see, that the total number of ks (massive particles), particles and holes (massless 
excitations) is always even. (It has to be noted, that the above mentioned freedom in 
choosing Jp introduces an ambiguity in the description of certain states: depending on the 
choice of Jp, in the same state the number of massless excitations can be different. This 
ambiguity drops out, however, of physical quantities, like the energy and momentum which 
do not depend on the choice of Jp.) The particles and holes being defined, it is convenient 
to write z(A) in the form 

2(A) = — | (sin _1 (A„ - iu) + sm~\A + iu 
2ir [ ^ 

I n(fc) j 

^E^i(A- sin kj) - - J2 V2 (A - A f ) 

- ^ E ^(A - A p ) + 1 £ ^ 2 (A - A*) J . (3.22) 
The density of the As is given by the derivative of the z: 
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i.e. 



«<A) = ^, (3.23, 



1 n ( k ) -1 

cr t (A) = ctq(A) - — 53 *i(A - sin k,) - — 2 ^ 2 (A - A/) 

1 53^ 2 (A-A p ) + -J^53^ 2 (A-A,). (3.24) 



2ttN^ y 2irN , 

P h 

Calculating the sum over the Fermi sea by means of Q3.4Q we arrive at 
a t (A) = t7 (A) - 53 Ki(A - sinkj) 

71 3=1 

1 E^(A-A p ) + -i-53K 2 (A-A /l ) 



2ttA^V 2vriV 



:^(A - 5 + ) + ——— 7 ^ K' 2 (A - B- 



27r2AN 2 a t (B+) zx ' 2iT2AN 2 a t (B- 

B+ 

- J K 2 (A-A')a t (A'). (3.25) 



B~ 

where the limits B + and B~ defined by the equations 



z{B ± ) = ±^ , (3.26a) 
which are equivalent to 

J „( A ) = / .«(A) = i {l - „- "W + ^-W + M j . (3 . 26b) 

As ( |3.25 ) is a linear equation, the density is a sum: 

<r t (A) = a b (A) + -^a w (A) + ^a {p} (A) + ^a {h} (A) + ^ /sc (A) , (3.27) 

where all terms (the contributions of the bulk, massive particles, the particles and holes, 
and the finite size corrections respectively) satisfy equations of the type 

B+ 

x{A) = 1 X (A) - — f K 2 (A - A')x(A') , (3.28) 
Ztt J 

B- 

where the inhomogeneous part, X X (A), can be written as: 
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a b (A): J 6 (A) = (7 Q (A) , (3.29a) 

a {k} (A) : J {fc} (A) = -— £ 2^ (A - sin%) , (3.29b) 

^ j=i 

<x w (A): J {p} (A) = -i-^ir 2 (A-A p ), (3.29c) 

<t w (A) : X {M (A) = +i- £ AT 2 (A - A,) , (3.29d) 

ff/ «(A) : I /<C (A) = -s^AWA - B + ) + iS^gzy *3(A - IT) . (3.29e) 

The equations for the holes and particles can be given in terms of the density: noticing, that 
the particles and holes will condense in the vicinity of the Fermi points, based on ( [3.21 ), 
( p.26| ) and ( |3.23| ) we may write: 

27T 



2vr / a t (A) = ^AJ+, (AJ+, = J+ h -J F = half - integer) , 

B+ 

2vr I «r t (A) = ^J- /h (AJ; A = -J F - J p/h = half - integer) . (3.30) 



p/h 



Here pj h is either p or /i, and the upper index + or — refers to the side of Fermi sea to which 
the particle or hole is near, and the 2ix factor is introduced for later convenience. Equations 
of the type ( |3.28|) in general (unlike the half-filled-band case) can not be solved in a closed 
form, neither the Eqs. ( |3.30| ) can be given in a more explicit way, nevertheless they can be 
handled in the u — > limit. 

4- Unbound electrons 



The unbound electrons are described by Eqs. (|2.4|) from which it is possible to eliminate 
the sum over the As using the summation formula Q3.4| ) leading to 

n(A) 

2nlj = Nkj — ^ (fi(sm kj — \ a ) — 

a=l 

- V?i(sin kj - A p ) + <Pi(sin kj - A h ) - 

V h 

' -K 1 (A-B+) + -— 1 K^A-B-)- 



2ANa t (B+) y ' 2ANa t (B-) 
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N / ^(A-A> t (A'). (3.31) 



The integral in ( |3.31|) can be transformed in the following way: if a function satisfies ( p.28|) , 



than it satisfies also the integral relation 
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B+ 



( B~ 



<^ m (£-A):r(A) 



B- 



+ J I <^ m (£-A)x(A) + 

y-oo b+ 

oo B+ 

+ J ^m(C- A)X X (A) - | ^ m+2 (e-A)x(A), 



(3.32) 



as it can be checked by calculating the convolution of ip m with x of (|3.28|) . It is easy to show 
by iteration, that this relation is equivalent to 



B+ 



A) z(A) 

(B- 



+ / 2 tan -1 tanh 



\-oo B+ 



TT(g- A), 

4w 



x(A) + / 2 tan -1 tanh 



^ - A) 
4u 



Xx(A) • (3.33) 



After substituting ( |3.33 ) into ( |3.31| ) and evaluating explicitly some of the integrals of the 
1(A) one arrives at 



i r i 7r(sin fc^ — A) 
27rJj = N Ikj- / 2 tan" 1 tanh -A -i -a (A) 



+ AT 



-oo 
oo ' 



J + J \ 2 tan -1 tanh 

\-oo £+ 



4m 

7r(sin kj — A) 



4m 



**(A) 



^2 tan' 1 5 J - + 

Q = l jV = 1 



^ i , 7r(sin fc 9 - — A D " 

J2 2 tan" 1 tanh - v 



7T 



4w 
1 



/ sin kj — sin kj> \ 

A ^ / 

+ 2 tan -1 tanh 
h 

1 



+ 



7T 



7r(sin kj — A/-, 
Au 
1 



2ANa t (B+) 2u cos h *(™^-* + ) 2 ANa t (B-) 2u cos h <!i^zgl) 



(3.34a) 

(3.34b) 

(3.34c) 
(3.34d) 
(3.34e) 



2h 



For a general it we find this is the simplest form of the equations for the unbound electrons, 
but these equations further simplify in the scaling limit. 

5. The energy and momentum 

To calculate the energy, in Q2.7|) we evaluate the sum over the A by using again the 
Euler-Maclaurin formula 



E = -Nut - (2t(cos kj -u)-fi) + J2 ^o(A p ) - £ e (A h ) 



2ANa t (B~ 



(3.35) 



Iterating a t in the integral by means of (|3.24 ) one arrives at 
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B+ 



E = -Nut + N J e(A)<7 (A) + 

B- 

+ V -(2t(coskj -u)-fi)- — I e(A)Kx(A- sin fa) 

2lX J 



J \ B- 
P h 

\ ^'(B + ) + l - -e'(B-) , (3.36) 

where e(A) satisfies Eq. (|3.28| ) with the inhomogeneous part: 

e(A): I e (A)=e (A). (3.37) 

(Note, that e ^ e r as the limits in the integral equations ( |3.28p and ( |3.12D are different.) 
Next, to calculate the integral of e x K\ we substitute the derivative of (|3.33| ), and obtain 

B+ 

E = -Nut + N [ e (A)<r (A) - (3.38a) 



B~ 

1 



\2ANa t[ B-)^ ~ MNMB-f'^ ) + ^ 
+ i>(A p )-$XA h ) + (3.38c) 

P h 

+ S ( - (2t(eos k s -«)-„)- ± J I „L,. A) £ p(A) ) + (3.38d) 
+—(/+/] — r z AT g(A). (3.38e) 

\-oo B+ / 2m 

In this expression the first and second terms contain the bulk contributions, the third term 
is connected with the finite size corrections (but also the second term contains finite size 
corrections due to the deviation of from ±B), and the rest is the excitation energy. 
The momentum is given by (|2.8|), which after summation over the Fermi sea reads: 



p = % E^-E^ + E^-E^ • (3-39) 

\ 3 a P h J 

As in (|3.30|) the particle and hole rapidities are connected to = ±J p /h — Jf type 

expressions, and as Jp may have more values, we write 

j a P h JV 

-E + E + - + -jf( An - An ~) . ( 3 - 4 °) 
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with 



An + = n(p + ) — n(h + ) , An = n(p ) — n(h )) . (3-41) 

Here the X^ ± ' ) means summations over the particles and holes near ± Jp, with n(p ± ) resp. 
n(/i ± ) being the numbers of the particles resp. holes near the Fermi points ±Jp. As the last 
term can give an infinite contribution in the scaling limit (where the momentum must be 
devided by a), we have to redefine the lattice so, that this term be equivalent to zero. This 
can be done, if the bandfilling n = 2N r /N is a rational number. Suppose, that n = l/q, 
where I and q are relative prime numbers, and that v and r\ are the smallest integers satifying 
v = Aqn/l. If in the redefined lattice v lattice sites form one elementary cell, than the point 
T[N T /N = nn/2 of the original Brillouin zone will be transformed into the origin of the new 
Brillouin zone, so the last term of ( 3.40|) can be dropped. (The redefinition of the lattice 



has consequences on the parities of some numbers: as the lattice must consist of an integer 
number of elementary cells, iV and N r must be integer multiples of 4 and 2, respectively) 



IV. EQUATIONS OF THE EXCITATIONS IN THE SCALING LIMIT 



1. The massless excitations 



In Ref. \'\ 1] the spectrum and structure of the massive excitations have been investigated. 
It has been found, that in order to keep the mass gap (m ) finite, relations (1.2) have to be 
obeyed. Now we examine the behaviour of the excitations connected with the A distribution 
in the limit fling ). 

First consider ( |3.30| ). As in the limit ( pL.2| ) N — ► oo, we divide these equations by a: 



- I *M 

a J 

B+ 



2ttAJ 



p/h 



-7 

a J 



2nAJ, 



p/h 



(4.1) 



P/h 



As the r.h.s. is finite, in order to have the l.h.s. finite too, 

^/h-B ± = 0(a). 

According to ( [L.2[ ) 



7rm 



u sin(7rn/2) exp 



7rsin(7rn/2) 
2u 



(4.2) 



(4.3) 



As the solutions of ( |3.28| ) change on the scale of u, the integrands in ( |3.30|) are constants on 
a scale oc a, and we may write: 



2TT<J t (B 



2ttAJ 



p/h 



2wo- t (B- 



B~ — A 



p/h 



2ttAJ 



p/h 



a L a L 

Finally, as a^B^) = a r (B) + 0(A _1 ), in the scaling limit, one obtains: 



2tt(jJB) 



A 



p/h 



B- 



2ttAJ 



p/h 



L 



2ncr r (B) 



B~ - A 



p/h 



2ttAJ, 



p/h 



L 



(4.4) 



(4.5) 
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These equations can be considered as the secular equations of the massless excitations in the 
scaling limit. They look very similar to the secular equations of free particles, but they are 
not exactly the same, as B + and B~, i.e. the location of the Fermi surface in the rapidity 
space depend on the number of excitations. 

2. Massive particles 

Now we consider the scaling limit of ( |3.34j ). The first term on the r.h.s. (line ( 3.34a|) ) 
can be evaluated to give 



i n -1 i 7f(sin kj — A) ,.. 4 _ f 7r 

kn — I 2tan tanh cr (A) = — vwexp \ 

Au n I 2u 



(4.6) 



This, in the limit (|1.2|) disappears for any n < 1. 
In the second term (line ( |3.34b| ) ) 



B~ oo\ 

j • j 2tan~ 1 tanh 

-oo b+ J 

a 



7r(sin kj — A) 



4?i 



\ -oo B 



2 exp 



*(A) = 
7r(sin kj — A) 



2u 



crt(A) . 



(4.7) 



As in the limit (|1.2j ) B ± = ±B + 0(1/N) (this can bee seen comparing ( |3.8b|) and ( |3.26b| )) 
and a t = cr r + 0(1/N), we may replace B ± by ±B, resp. a t by o r (the errors introduced this 
way are proportinal to a/u resp. a). Thus, after some simple manipulations we have 



( B- 



N 



i , 1 , 7r(sin kj — A) ... 

+ / | 2 tan" 1 tanh — / -a t (A) 

Au 



\-oo B+ 



L ~a 6XP {~lf} / 6XP {^} ^ + A) Sinh 



77 sin kj 
2u 



(4.8) 



Finally using that ( [0, glf 

7rA 



/r, * \ 2u [W Au 2 Hi 

exp <; -— f> <r r (5 + A) = —J-a(B) + 



2u 



7T v e 



7T^ 



(4.9) 



and 



lim«r(S) = ^lima (5) 



(4.10) 



to leading order, after substituting the value of B ( 3.14Q , we find that the contribution of 
this term to the r.h.s. of (12.4) is 



f 7T sm kj 
Lm sinh < 

2u 



(4.11) 
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with m given by (|1.2|) . 

The next two terms (line ( 3.34cj ), the contribution of the As and ks) survives the limit. 



The contribution of the holes and particles in the A distribution (line ( |3.34d| ) tends to a 
constant 

^(An+ - An~) , (4.12) 

while the last two terms (the finite size effects, line ( |3.34e|) ) disappear. 

In the parity prescription for Ij there is a parameter, n(A), tending to infinity in the 
scaling limit. It can be replaced by finite numbers in the following way: first of all, as 
n(A) = rip + n(p) — n(h) with rip being the number of elements in the Fermi sea, and that 
the total number of 'excitations' n(k) + n(p) + n(h) =even, n(A) = rip + n(k) (mod 2). It is 
to be noted, that n F = N r + 5, i.e. it can be both even and odd, but due to the redefinition 
of the lattice N r is always even, thus we have 

/ _ n(A)+n(fe) + V odi) . (4.13) 
Collecting all the terms, and introducing the notations 

^ = !L lr 1 ' ("(«)=*(*))» ^ = ("(*) ="(*))> (4-14) 

we conclude, that the massive particles are described by the equations 



3 = 4 r 



Lp( Kj ) = 2nl> - £ {^-^-) + E Stan" 1 (^^^J , (4.15a) 
An + — ArT 



j 



2n(x) + 2n(«) + 25- {An + - An 
A 

nin) — 2n(x) + 2m n 

+ — (mod 1 



(mod 1~! 



4 2 



E^an- 1 (^-^j =27rJ a + X:2tan- 1 (^^J , (4.15b) 

modi) ) , 




n{n) + n(x) + 1 
V a = 2 

with 

p(n) = m sinh k (4.16) 



Eqs. Q4.15D are the secular equations of the massive particles. These equations have the same 



structure as the corresponding equations of the half-filled band case, with one significant 
difference: the parity prescription for the quantum numbers 7j depend not only on the 
number of excitations, but also on the parameter 5. (For later purposes in the last row of 
( |4.15a| ) we give the Ij in terms of other quantum numbers n and m defined in ( |6.11 ).) 
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V. ENERGY AND MOMENTUM IN THE SCALING LIMIT 



1. The energy 

Now we calculate the scaling limit of G3.38|). First we calculate the contribution of the 



unbound electrons (lines (|3.38d|) and ( |3.38e|) ). The expression in ( |3.38d|) is exactly the same, 
as the energy of a massive particle in the half-filled band: 



1 °° 

— (2t(cos k — u) — fi) / 

277 J 



7T 



2u cosh " (si ° fc ~ A) 

2u 



— \/uexp{ X chn . 

ix I 2u J 



(5.1) 



These terms disappear in the scaling limit. In ( |3.38ej) we may replace B ± by ±B and e(A) 
by e r (A), and we arrive at 




7T 



2u C osh " (sin 9 fc ~ A) 

2u 



s(A) 



u 



exp 



txB 
~2u 



exp 



7rA 

2u ) 



\ 



e r (B + A) cosh 



77 sin k 
2u 



(5.2) 



The integral on the r.h.s. can be evaluated ( |T4|, JTT 

4m 2 pif 



exp 



ttAI 
'~2u) 



s r (B + A) 



lime'JB). 

7T 2 V 2e«^o 01 ; 



(5.3) 



(The structure of this relation is the same as that of (|4.9|), the difference is due to ( |3.13| ).) 
Finally we find that the contribution of the massive particles is simply given by 



•31 1 



(5.4) 



where 



e(/c) = mo cosh k . 



(5.5) 



Next we consider the energy contribution of the particles and holes in the A distribution 
( gggg ). Due to (p!3Q 

e(A) =£ r (A)+o(5 ± T fi), (5.6) 

thus e can be replaced by e r . Moreover, as the e r changes on a scale oc it, but A^ 1 =|= £> oc a, 
we may linearize e r around ±B. Finally we arrive at 

Ee(Ap)-£^) = 

= E (+) 4(£)(a+ - B + )-£ (4 V r (a)(A+ - 5 + ) - 

V h 



-E 4(^)(A; " B-)+J2 e' r (B)(K ~ B-) + 

P h 

+£' r (B)(An + (B + -B) + An-(-B~ - B)) . 



(5.7) 
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This way the contributions of the particles and holes is expressed by the quantities known 
from fl45|) : 



2nv fv^(+) 



(e (+) a 4- E (+) aa + + E M a j p - - E M A4- } + 

\ P h V h ) 



/ , a r (B)(B+ -B) A a r (B)(-B- - B)' 
+2irv An + rK JK '- + An~ rK A ^ 



(5.8) 



where 



(5.9) 

We note that due to (|3.13|) , s' r (A) satisfies an equation of the type (|3.12| ), with e (A) replaced 
by £q(A). A consequence of this is that analogously to ( |4.1U| ) 

lim4(S) = -= \im£' (B), (5.10) 



which leads to 



v = l (5.11) 

in the scaling limit. 

The terms in line ( |3.38b| ), after replacing s'^B^) and a t (B ) by s' r {B) and a r (B), re- 
spectively, turn out to be 

~w (5 - 12) 

Finally we consider the terms in line ( |3.38a| ), where the integral is obviously a function 
of B~ and B + but it also depends implicitly on e(A). It can be expanded into a power series 
of (B + — B) and (B~ + B) and after a straightforward calculation we find, that 

B+ 



- Xirl + X I s(A)<r (A) = E r + ^2nL 



a r (B)(B+-B) 



a r (B)(B- + B) 



Collecting all the terms (|5.4|-|5.13|) we have 



7T 



E = E r + 

r 6L 



+ -2vrL 
2 



'a r (B)(B+-B) 



'a r (B)(B~ + BY 



+ 



( An+ * r (B)(B+-B) An _a,. m -B-- B y 



{e (+, aj p + - £ (+, aj+ + z ( ~>^ r -- E H Mr| + 

I v h v h ) 



2vr | v ^(+) 
L 



+ E e (^)- 



(5.13) 



(5.14) 
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Next one has to calculate the quantities a r (B)(B + — B)/a and a r (B)(—B~ — B)/a and 
to do this we use the conditions ( 3.8b ) and ( |3.26b|) . Taking the difference of these two 
equations, and keeping only the terms of order of 1/N, we arrive at the relation 

c Ub )( b+ - B ) + _ ar{B){B+ _ B) = 



2N 

where 



— (n(/c) + An + + An' + 25) , (5.15) 



C = Jp(A), (5.16) 

B 

with p determined by the equation 

1 1 f 

p(A) = ~^K 2 (A -B)--j K 2 (A - A')p(A') , (5.17) 

—B 

and at an analogous relation for {—B—B~). In the u —>■ limit Eq. ( |5.17| ) can be transformed 
into a Wiener-Hopf type equation and the integral (|5.16|) can be calculated: 

C=l-V2. (5.18) 

This way we find, that 

a r (B)(B+ -B) n(p + )-n(h + ) n(«) + 3An+ + An" + 25 
a L 2y/2L 

Substituting this and the analogous expression for (— B — B~) into (|5.14j ) yields 
2tt 1 f (nU) + 2An+ + 2/An- + 25\ 2 f An + - An" \ 2 



+ T2U^ — 2 ; + 

4{e ,+ v-e ,+, aa + -^}+ 

4{i: h ^--i: w a4--^}. (5.20) 



Note, that the values of the expressions in the last two curly brackets are always nonnegative 
integers. 

2. The momentum 

The momentum on a lattice is a dimensionless number, while in a continuum it has the 
dimension 1 /length. In the continuum limit the momentum, P, is obtained by taking the 
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limit of Pi j a (with P% being the lattice momentum). Now P\ is given by ( |3.40|) which after 
dividing by a and substituting Eqs. Q4.15|) yields: 



that is 



+E 

V 

-E 



(+) 2ttAJJ 



L 



E 



(+) 2vrAJ, H 



(-) 2ttA J" , ^(-) 27rAJ h - 



7T 



H — p(n(/«) + 28) (An + - An" 



(5.21) 



p = EMs- 



2tt /n(/t) + 2An+ + 2An- + 28\ ( An + - An" 
T [ 2 J I 2 



+ 



+ ^(e ,+, a^-E <+, a^-^ 

(An-) 2 



L 

2n 



-9^E h aj p --E^aj ; 



(-) 



(5.22) 



VI. INTERPRETATION 



The secular equations, Eqs. (|4.5|) (for the massless sector) and Q4.15D (for the massive 



sector) do not admit an immediate interpretation in terms of scattering states of a massive 
and a massless SU(2) doublet as this has been the case for the scaling limit of the HF 
Hubbard chain fIC| . We present below our interpretation of the results obtained in the 
previous sections. 

(i) The ground state (the physical vacuum) is the lowest energy state of N r bound pairs of 
bare particles. This state can be excited by intoducing 'particles' and 'holes' into the 
ground-state distribution of the bound pairs and/or introducing unbound electrons. 
The former set of excitations resembles a set of massless particles, while the unbound 
electrons are massive. The two sectors do not completely decouple: 

- the total number of (massive and massless) excitations must be even, this is a 
consequence of the parity prescriptions the quantum numbers must obey, and the 
requirement, that the Fermi-sea is symmetric; 

- the parity prescription for the quantum numbers of the massive particles depends 
on the state of the massless sector; 

- in the energy the numbers of massive and massless particles are nonlinearly coupled. 
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(ii) As we have noted, the definition of the Fermi sea we used is only one of several (actually 
infinite) possibilities: the same state can be described in terms of a larger or smaller 
Fermi sea. We may choose 5 even to be a function of of the described state, but the 
rule 

ti(k) +n(A) = n r (A) +5 (mod 2) (6.1) 

must be obeyed. It is obvious, that n(p) and n(h) (so An + and An - ) depend on the 
definition of the Fermi sea, but the quantity entering the energy and the momentum 

ti(k) + 2An + + 2An- + 25 = N e - 2N r (6.2) 

is independent of the definition of 5, as it gives the deviation of the actual bare particle 
number from that of the reference state. It is easy to see, that the other quantity which 
enters into the energy and momentum 

An + - An' (6.3) 

is also independent of 5, hence no physical (measurable) quantity depends on our 
choice of the Fermi sea. It is, however, unavoidable to make a definite choice in order 
to be able to define the massless particles. 

A simple consequence of the above arbitrariness is that we can formally 'minimize' the 
coupling between the massive and massless sectors: we can choose a 5 (obeying ( |6 . 1| ) ) 
such that 5' = u(k)/2 + 5 have the smallest modulus. There are four inequivalent 
values for 5': or 1, if n(/c) is even, and +1/2 or -1/2, if is odd. We note, that 
Eq. (|6.2|) has now the form 

2An +l + 2An-' + 25' = N e - 2N r . (6.4) 

The 'minimizing' of the coupling between the massive and massless sectors corresponds 
to choosing the Fermi level with 2An +l + 2Arr' as close to N e — 2N r as possible. It 
is worth to mention that in the HF band case 2N r — N e = N — N e is measured by the 
number of massless particles only, in this sense the above choice of 5' mimicks most 
the HF case. 

(iii) To specify a state it is not sufficient to give the number of particles (and the quantum 
numbers) in the two sectors, but one also needs the value of 5 or 5'. It is not hard 
to see, that any state with given n(ft), n(p) and n(h) can have both 5 = and 5=1 
(if we use the convention ( |3.18| )), or equivalently can have both 5' = and 5' = 1 for 
ti(k) even, and 5' = ±1/2 for ti(k) odd. In this sense 5 (or 5') is a free parameter. 
States differing only in the value of 5 (5') are not degenerate, as 5 appears in the parity 
prescription for Jj in (|4.15a|) , and it appears explicitly in the energy (|5.20|) too. 



(iv) The massive sector described by Eqs. ( [4.15|) consists of (relativistic) particles in the 
doublet representation of SU(2). Their contribution to the energy and momentum 
is not simply the sum of the individual contributions but there are terms both in 
the energy and momentum, which depend on the number of massive particles in a 
nontrivial way. (As these terms contain data on the massless sector too, we interpret 
them coupling between the two sectors.) 
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Any solution of Eqs. ( |4.15| ) corresponds to a highest weight state of SU(2) with S 2 = 
1(1 + 1), S z = m, I = m = (ti(k) — 2n(x))/2. All the other members of the / = 
(ti(k) — 2n(x))/2 multiplets are degenerate with this state, and can be obtained by 
the action of the a~ = J2i c /i c i,T operator. 



The BA Eqs. ( 4.1 5| ) are of the familiar type and this makes possible extract the two 



particle scattering matrices in the usual way |17-19[]. For this we have to choose 



those solutions of Eqs. (|4.15|) which correspond to the triplet and singlet states with 



an empty massless sector. To avoid ambiguities due to the freedom in choosing the 
Fermi level one should first define the state with no massless particles. It is natural to 
consider the massless sector empty, if its energy contribution is zero. This corresponds 
to m = n = (see the next point), and then the scattering matrix of the massive 
doublet is given by 

S Ak = - exp U — — — , (6.5) 

I V 7T / J — ITX 

where I resp. fl are the identity resp. permutation operators acting on the spins of the 
two particles: 

ri a 2 — a s n CT i ,7 2 = 5 5 (6 6) 



(v) The massless sector described by Eqs. ( |4.5|) . looks very much like as if it consisted of 
free particles. This picture, however, cannot be taken at face value as the definition 
of the particles (and holes) is not unique. Also as both the energy and momentum 
contain the collective terms it looks as if there were interactions. One cannot, however, 
extract phaseshifts from Eqs. 



The energy and momentum can be described in terms of a CFT. If the massive sector 
is empty the energy momentum dispersion shows a tower structure: the central charge 
is c = 1, and the apexes of the towers are located at 

6L L L 



with 



x 



n,m 



2 

where the integers n and m are 



- (n 2 + m 2 ) , s n , m = nm , (6.8) 



n 



t » j- <-\ ( An + — An \ 

(An + + An' + 5) , m = I . (6.9) 



This leads to the conformal weights (the notations are those of [EH]) 
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A = ^(n + m) 2 , A = i(n-m) 2 . (6.10) 

(We have to note, that even if the massive sector is not empty, the contribution of the 
massless sector is of the tower structure, just with 

(nU) + 2An+ + 2An~ + 2S\ ( An+ - An - \ 
n= { 2 J' m= ( 2 )■ (6 - U) 

In this case, however, due to the coupling of the massive sector, the apexes of the 
towers can not be interpreted as anomalous dimensions and spins. For n(/c) ^0n and 
m can be half-integers too: 2n,2m = n(n) (modi).) The conformal weights fl6.10|) 
correspond to a SU(2)xSU(2) symmetric CFT. It should be interesting to directly 
identify the corresponding SU(2) xSU(2) multiplet structure. 

(vi) The spectrum of the limiting model is the same both in the NHF and the HF band 
case. Actually the complete energy and momentum are 

„ _ / % n 2tt ( n 2 + m 2 \ 2tt / , \ 

E = E r + £ efo) - — + — [-^—) + Y ( + ' (6 - 12) 

p = J^PiKj) + ^nm + ^ (u + - u~) . (6.13) 

J 

The towers are obviously the same as in the HF case, and so is the contribution of the 
massive sector, as the parity prescription for the quantum numbers in (|4.15|) (i.e. the 
quantization of the momenta of the massive particles) when expressed in terms of n 
and m agrees with that of the HF case. We note, however, that this way we have only 
shown that the points of the two spectra coincide, but we can not say anything about 
the degeneracies of the single points. 

Finally let us briefly present the results for the ground state energy of a finite density 
state (i.e. when the number of excitations is macroscopic) of the limiting model. As 5 can 
be choosen to minimize the coupling between the sectors, one can discuss the massive and 
massless sectors separately. (The coupling due to the parity-prescription for Jj and the 
explicit presence of 5' in the energy of the massless sector may actually give a contribution 
0(1/L) in the energy density.) Since then the massive sector is described by the same 
equations and energy-momentum dispersion as in the HF-band case, one can literally take 
over the results obtained there |T(J and we do not reproduce the formulae here. 

The treatment of the massless sector is extremely simple. After choosing the 5' as 
described in the previous section, the energy density associated with the massless sector is 
just given as 

£(q)=7t 6 2 , (6.14) 

where g is the density of massless particles. Eq. (|6.14j) coincides with the energy density of 
noninteracting particles. 
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The total free energy density (defined as the Legendre transformation of E(q) with a 
chemical potential v) will be simply the sum of the massive and massless contributions just 
as in the HF case. The decoupling of the two sectors depends apparently on our choice of the 
Fermi level. We can understand this as follows. There are two energy contributions controled 
by two independent parameters. One is the energy of the massive particles, related to their 
relativistic dispersion, this is determined by n(/c). The other contribution is a (quadratic) 
function of iV e — 2N r , i.e. the deviation of the actual particle number from that of the 
reference state. If we make the Legendre transformation in this two variables, we shall get 
two independent terms. If we choose some combinations of n(n) and N e —2N r as independent 
variables, the free energy will have terms depending on both chemical potentials. (In the 
HF-band case 2N r — N e = N — N e is measured by the number of massless particles, thus 
the separation of the two sectors is obvious. For the NHF-band case measuring N e — 2N r 
can be done in different ways. After having choosen the 5', N e — 2N r ~ 2Lg. The factor 
2 in this equation appears as g is the density of bound pairs and this is the reason for the 
discrepancy of (|6.14| ) from the analogous expression in the HF case.) 

In conlusion the above result provide strong (although indirect) evidence that the scaling 
limit of the NHF Hubbard chain is the SO (4) symmetric CGN model. 

VII. DIFFERENCES BETWEEN THE HF AND NHF CASES 

In the formulae ( |1.2D one can take the n — > 1 limit suggesting that there is a smooth 
limit from the NHF to the HF case. In this section we list some differences between the HF 
and NHF cases, which clearly show, that in spite of the 'continuity' of ( |1.2| ) in n the scaling 
limit of the NHF and that of the HF case can not be related in a trivial way (in other words 
the n — > 1 and the scaling limit do not commute in a obvious way). 

• In the HF case the ground state is well defined, and in the description of the bound 
pairs (A) the Fermi surface plays no role: the Fermi points in the rapidity space (5 ± ) 
are at infinity. As a consequence, there are no particles, only holes connected with the 
A distribution. These elementary excitations (we call them massless particles) carry 
an SU(2) isospin and the states are isospin eigenstates often characterized by a set 
of complex variables. In the NHF case the Fermi surface is at some finite rapidity, 
there are both particles and holes, but in the limiting theory there are no states with 
complex As, as those require an infinite excitation energy. 

• In the HF case the two kinds of excitations have two different internal degrees of 
freedom (spin resp. isospin), and the ground state is a singlet of both. In the NHF 
case one of the excitations carries the spin, but the isospin of the state cannot be 
uniquely attached to one or the other kind of excitation. The reference state is a 
spin singlet, but its isospin (N — 2N r )/2 —>■ oo in the scaling limit. (The difference 
between the isospin of an excited state and that of the reference state is, however, 
finite: A/ 3 = -(n(«) + 2An + + 2An~ + 25) /2) 

• In the HF chain the massless particles are uniquely defined, in the NHF case - 
although the state can be given uniquely — the definition of the particles and holes 
connected with the A distribution (massless particles) is not unique. 
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• In both cases the conformal dimensions of the massless sector correspond to an 
SU(2)xSU(2) symmetric CFT. In the case of the HF chain an SU(2) symmetry corre- 
sponding to the gapless excitations is already present in the lattice model, and this with 
the separation of the left and right sectors develops into an SU(2)xSU(2) symmetry. 
For the NHF chain in the lattice model there is no such SU(2) symmetry, which could 
be identified in the massless sector and the SU(2) xSU(2) symmetry (whose existence 
is indicated by the conformal weights) appears only in the scaling limit. 

• While for the HF case the conformal weights A and A are connected to the right resp. 
left going massless particles, for the NHF case both A and A depend on both the right 
(+) and left (— ) part of the massless sector. 



APPENDIX A: 



We sketch here briefly the 'naive' continuum limit of the less than half filled Hubbard 



chain. A more detailed description of the procedure is given in Ref. pi 
First we write the (II . If) Hamiltonian in the form 



N jj N 

H = ~tJ2 E ( C t* C i+l,v + h - c ) + 17 E ("i,T + h hif 
i=la=t,l ' / i=l 

A NU 

i=i 4 

with h = fi — U. 

We have seen earlier, that in order to avoid divergent terms in the momentum we have to 
redefine the lattice so, that in the new lattice v sites form one elementary cell, {y together 
with an other number 77 is defined so, that if the bandfilling is n = l/q, where / and q are 
relative prime numbers, than v and rj are the smallest integers satifying v = Aqrj/l.) We 
define the operators 

1 p 2n 
(n) = ^=Y, e ~ lkaJc »(n-D+^ a = l,2,...,u, a^t,!, k = — . (A2) 



In terms of these the original Fermion operators are 

c u(n - 1)+j , a = Va^ e ikaj <p a ^(n) j = 1, 2, . . . , v. (A3) 

a=l 

The part of the Hamiltonian quadratic in the field operators now reads 

N' v 

EEE pa ( h ~ ^ cos ka)(f)^ a (n)(f) a>a (n) 

n=l a a=l 

N' v 

"«* EE E <M{Mn + l)-Mn))e lkp 

n=l cr a,/3=l 
N' v 

+«* E E E <M (Mn) - faAn ~ 1)) e~ ika , (A4) 

n=l cr a, /3=1 
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with N' = N/ v being the number of new elementary cells. This expression shows that the 
procedure leads to a meaningful result only, if for some a = «o and a = v — «o 

h - 2t cos ka = 0, (A5) 

otherwise all the new fields become infinitely massive if t — > oo. Obviously 

—J- = M= -y) ( A 6) 

what implies, that 

a = V , ( A7 ) 

with r; defined above. It is clear, that all the modes a ^ rj,i/ — rj can be omitted, as they 
become infinitely massive in the continuum limit. 
The length of the chain is 



L = Na = vaN' (N' = int.) 
and the continuum limit is defined as 

a — ■> 0, N' — > oo, L = fixed 

with the continuous variable 

x = vain — 1/2), dx = va. 



(AS) 



In this limit 



N' 

E 



dx 



va 



and 5 ntn f — > va5(x — x). 



If we introduce 



then 



<Pr h a{n) = iPi,<t(x), 4>u-r,A n ) = VvO*)' 



(A9) 



(A10) 



(All) 



(A12) 



Applying (|A3|) , (|A11| ) and ( |A10|) to the Hubbard Hamiltonian ( |Al"| ) one obtains in the naive 
continuum limit the following Hamiltonian density: 



H(x) = (2at sinTm/2) - £ ^ 5 d x ^ a + 



a=l 



+ U 



era 1 V a / 



where 
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i 
-i 



(A13) 



(A14) 



(A15) 



and u = f//4tsin7rn/2. (Note that this differs from the definition used in the bulk of the 
paper as there u = \U\/4t.) 
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